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ABSTRACT 

A  review  of  procedures  for  designing  acceptance 
sampling  plans  and  process  control  schemes  according 
to  economic  criteria  is  given.  The  use  of  nonsta- 
ttstlcal  criteria  to  design  these  procedures  has 
bean  the  subject  of  considerable  research,  and  a 
number  of  different  process  models  and  procedures 
have  been  developed.  A  comparison  of  the  major  ap¬ 
proaches  to  model  formulation  and  design,  and  a  dis¬ 
cussion  of  the  practical  implementation  of  these 
techniques  is  the  primary  focus  of  this  paper. 

INTRODUCTION 

Acceptance  sampling  and  process  control  are  the 
primary  statistical  techniques  used  In  a  quality 
assurahee  program.  These  techniques  have  found  wide 
application  in  industry,  particularly  in  manufactur¬ 
ing,  In  the  following  areas:  incoming  material 
Inspection,  surveillance  of  the  production  process* 
estimation  of  lot  or  process  characteristics,  process 
capability  analysis,  and  finished  product  quality 
auditing. 

Traditionally,  acceptance  sampling  and  process 
Control  schemes  are  designed  with  respect  to  statis¬ 
tical  criteria.  For  example*  one  may  choose  the  sam¬ 
ple  else  and  acceptance  number  for  e  single- sampling 
plan  for  attrlbutea  so  that  tha  operating  character¬ 
istic  curve  passes  through  (or  near)  two  points 
specifically  selected  to  give  certain  probabilities 
of  lot  rejection  at  specified  levels  of  lot  or  pro¬ 
cess  quality.  Similarly,  one  may  design  a  control 
chart  so  that  the  power  of  the  chart  to  detect  a 
particular  shift  in  process  quality  and  the  probabi¬ 
lity  of  false  alarms  ere  equal  to  specified  values. 
While  the  traditional  approach  often  produces  accep¬ 
table  results,  It  is  also  possible  to  design  accep¬ 
tance  sampling  and  process  control  schemes  with  re¬ 
spect  to  economic  criteria.  This  has  considerable 
Intuit iva  appeal  since  there  has  been  an  increasing 
emphasis  on  quality  costs  in  recent  years.  Further¬ 
more,  the  use  of  these  techniques  have  direct  econo¬ 
mic  consequences  In  tbet  one  Is  balancing  the  costs 
associated  with  sampling,  testing  and  process  sur¬ 
veillance  against  internal  and  external  failure 
coats.  Sines  tha  daslgn  of  the  procedure  effects 
these  costs,  It  la  logical  to  consider  this  design 
from  an  aconomlc  viewpoint. 


During  the  last  20  years,  considerable  re¬ 
search  has  been  devoted  to  the  use  of  economic  cri¬ 
teria  to  design  acceptance  sampling  plans  and  pro¬ 
cess  control  schemes.  This  paper  gives  an  overview 
of  the  major  developments  and  approaches.  Some 
comparison  of  the  different  approaches  to  model 
formulation  and  a  discussion  of  the  practical  Imple¬ 
mentation  of  these  techniques  Is  also  given. 

ECONOMIC  MODELING  AND  PROCESS  CONTROL 

The  Shewhart  control  chart  is  probably  the  most 
widely  used  process  surveillance  device.  To  design 
a  control  chart,  we  must  choose  the  sample  size  (n), 
the  control  limit  (k,  the  multiple  of  o) ,  sod  the 
Interval  between  samples  (h  hours) •  Despite  the 
non-optimality  of  fixed  sample  size,  fixed  sampling 
Interval  procedures,  the  Shewhart  control  chart  has 
gained  widespread  use  because  of  Its  flexibility, 
simplicity  of  administration  and  the  additional  In¬ 
formation  about  process  performance  often  contained 
In  the  pattern  of  points  plotted  on  the  chart. 

Considerable  effort  has  been  devoted  to  deve¬ 
loping  economic  models  of  Shewhart  control  chart#. 
Theaa  models  usually  assume  that  the  process  is 
character iced  by  one  in-coat rol  state  that  represents 
the  mean  of  the  quality  characteristic  when  no 
assignable  causes  are  present,  and  s  *  1  out-of- 
control  states.  The  probability  model  chat  governs 
the  transitions  between  these  i  +  1  states  Is  called 
the  process  fsliurs  mechanism. 

•  1 

Costs  and  Measures  of  Effectiveness 

Three  cetegories  of  costs  are  usually  consid¬ 
ered  In  the  development  of  economic  models  for 
Shewhart  control  charts;  tbs  cost  of  sampling  end 
testing,  usually  of  the  form  a^  +  a2n,  the  costs  of 
investigating  and  possibly  correcting  action  signals, 
*and  the  coat  of  producing  defective  items  (Internal 
and  external  failure  costs) .  Thane  costs  arc  then 
combined  to  form  a  total  cost  per  unit  time  function. 
Tbs  general  approach  la  to  define  a  cycle  as  the 
length  of  time  T  during  which  the  process  begins 
operating  In  the  in-control  stats  and  eventually  re¬ 
turns  to  the  ln-control  stats  following  a  process 
adjustment. 
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If  C  la  the  coat  incurred  during  a  cycle,  then 
the  expected  cost  per  unit  time  is 

EU)  -  E(C)/E(T)  (i) 

Since  E(C)  and  E(T)  are  functions  of  the  design  par* 
enters  of  the  control  chart,  numerical  optimization 
(direct  search)  techniques  car  be  employed  to  find 
the  parameter  values  that  minimize  the  expected  cost 
per  unit  time. 

Basic  Process  Models 

The  fundamental  process  model  often  used  in 
economic  design  of  the  Shewhart  control  c..art  is  due 
to  (hmcan  ( 7 1 *  His  work  was  the  first  to  deal  with 
a  fully-economic  approach  to  control  chart  design 
and  to  explicitly  consider  the  optimization  problem. 
The  IXincan  model  assumes  a  single  assignable  cause 
(thus  a-1)  and  Chat  transitions  between  the  in- 
control  and  out-of-control  states  occur  according 
to  a  Poisson  process  with  intensity  X  occurrences 
per  unit  time.  Thus  the  length  of  time  the  process 
remains  in  control,  given  that  it  starts  in  that 
state.  Is  an  exponential  random  variable  with  mean 
1/A.  This  is  equivalent  to  a  process  failure  mech¬ 
anism  with  a  uniform  hazard  function. 

Tha  cycle  consists  of  four  periods:  (1)  the 
in-control  period,  (2)  the  out -of -control  period,  (3) 
tha  time  to  taka  a  sample  and  interpret  the  results 
and  (4)  tha  time  to  find  the  assignable  cause.  As¬ 
suming  that  the  process  continues  to  run  during 
searches,  the  expected  cycle  length  is 

K(T)  -  1/A  +  h/(l-0)  -  t  +  gn  +  D  (2) 

where  h  la  the  sampling  interval  (in  hours),  1  -  0 
Is  tbs  power  of  the  chart  for  detecting  a  specified 
out-of-control  stats,  n  is  the  sample  site,  g  is 
tbs  time  required  to  take  and  Interpret  e  sample,  D 
Is  the  expected  time  to  find  an  assignable  cause  and 

1  (J+l)h  (j+l)h 

T  m  f  «~At  X(t-Jh)dt  / f  «“Xt  Adt 
Jb  JH 

!  -  [1  -  (l+Ah)*‘Xh]  /  A(l-._Xh)  (3) 

,1a  the  expected  time  of  occurrence  of  the  shift, 
given  that  it  occurs  betweeen  the  and  (j+l)^ 

samples. 

Let  Vo  end  be  the  net  income  per  hour  of 
operation  in  the  in-control  and  out-of-control 
states,  respectively.  The  costs  of  investigating 
real  end  false  alarms  are  *3  and  el,  respectively. 
Since. there  are  E(T)/h  samples  par  cycle  and  oe'*h/ 
false  alarms  par  cycle,  on  the  average,  the 
expected  net  income  per  cycle  Is 

»<C)  -  V0(l/A)  +  (Sj+Sjn)  E(T)/h 
♦  V1[h/(1-S)  -  T  +  gn  ♦  D] 

-  »j  -  «3«."Xh  /  a-«*Xh)  <*) 


Dividing  (4)  by  (2)  gives  the  expected  net  in¬ 
come  pt-r  hut.  r ,  say 

E(I)  -  V0  -  E(L)  (5) 

where 

E(L)  -  Uj+a^/h 

+  Ja^[h/(l-0)  -  t  +  gn  ►  D] 

.  ,  *  -Ah  .  ,,  -Ah./ 

♦  *3  +  /  (1-e  )\ 

*  ll/x  +  h/(l-0)  -  T  +  gn  ♦  D]  (6) 

and  ■  Vq  -  is  the  penalty  coat  (par  hour)  re¬ 
sulting  from  production  in  the  out-of  control  state. 
Maximizing  E(I)  In  (5)  Is  equivalent  to  minimizing 
E(L)  in  (6).  There  are  three  decision  parameters, 
n,  k  and  h  (the  function  E(L)  depends  on  the  width 
of  the  control  limits  k  through  a  end  0);  end  nine 
user-supplied  constants  that  describe  the  production 
process,  a^,  S2»  S3,  S3,  *4,\,  g,  D,  and  $  tha  mag¬ 
nitude  of  the  process  shift  (which  Is  Involved  In 
determining  o  and  0). 

While  this  basic  process  model  could  be  used 
in  many  situations,  there  are  a  number  of  variations 
that  are  potentially  important.  The  first  of  these 
Is  the  Incorporation  of  more  than  one  out -of -control 
state.  Duncan  [8]  has  generalised  the  cast  model 
above  to  a  situation  where  there  are  ■  out -of -control 
states  (assignable  causes).  Knappenberger  and  Cran- 
dage  [13]  also  prssent  a  model  capable  of  treating 
several  out-of -control  states.  The  Duncan  model 
appears  somewhat  more  real let 1c  in  that  It  employe 
different  search  costs  for  different  assignable 
causes,  while  the  Kneppeaberger  and  Grandage  modal 
uses  an  average  or  composite  cost.  On  the  other 
hand,  thq  Knappenberger  end  Gramdege  modal  allow* 
continued  process  deter lore Cion  beyond  the  Initial 
shift,  while  the  Duncan  modal  allows  only  one  or  two 
shifts  to  occur.  Both  modsls  are  considerably  more 
sophisticated  in  terms  of  thalr  data  requirement a 
than  (6)  above.  Furthermore,  there  is  considerable 
evidence  that  a  single  assignable  causa  model  that 
matches  tha  true  multi-state  system  In  Important 
ways  will  be  s  satisfactory  approximation.  Specifi¬ 
cally,  if  we  let  the  single  out -of -control  state  be 
a  weighted  average  of  the  s  out-of-control  states 
with  the  weights  chosen  proport ionslly  to  the  proba¬ 
bility  of  occurrence  of  tha  state,  excellent  results 
ere  obtained.  In  practice,  there  is  probably  little 
reason  to  use  a  multiple  out-of -control  state  model 
in  most  situations. 

The  assumption  of  s  constant  hazard  function 
for  the  process  failure  mechanism  (an  exponential 
distribution  for  the  time  between  shifts)  is  criti¬ 
cal.  If  tha  occurrence  of  assignable  causes  can  be 
thought  of  as  random  "shocks"  disturbing  the  systma, 
then  the  probability  of  s  thlft  occurring  within  any 
smell  interval  of  time  is  proportional  to  the  length 
of  the  Interval,  and  the  exponential  model  is  prob¬ 
ably  appropriate.  However,  If  assignable  causes 
occur  as  s  result  of  the  cumulative  effects  of  vibra¬ 
tion  or  host.  Improper  set-up,  or  excessive  etressee 
during  start -up,  as  Is  tha  case  with  many  production 
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processes,  then  the  exponential  model  may  be  entirely 
Inappropriate.  Distr Ibutions  with  either  increasing 
or  decreasing  hazard  functions  may  be  more  reasonable. 
His specification  of  this  aspect  of  the  process  model 
may  have  very  serious  economic  Impact.  Unfortunately, 
It  Is  not  straightforward  to  incorporate  other  process 
failure  mechanisms  in  (6).  A  discrete-t ime,  single¬ 
cause  model  that  allows  free  choice  of  the  process 
failure  mechanism  has  been  developed  by  Baker  [2], 
but  as  it  assumes  that  a  sample  is  taken  at  the  end 
of  each  period*  the  optimum  inter sample  interval  can¬ 
not  be  explicitly  determined. 

The  assumptions  in  (6)  that  the  process  contin¬ 
ues  In  operation  during  the  search  for  ar  assignable 
cause  and  that  the  co«t  of  repair  is  not  Included  are 
unrealistic  for  many  processes.  Fortunately,  (6)  may 
be  easily  modified  to  consider  these  possibilities. 

If  Dq  and  D2  are  the  expected  search  times  for  false 
and  real  alarms,  respectively,  and  if  &  is  the  cost 
of  repair,  then  asmimlng  that  the  process  Is  shut 
down  during  the  search  leads  to  an  expected  net  in¬ 
come  per  unit  time  of 

E(I)  -  VQ  -  E(L)  (7) 

where 

I(L)  -  j(«1+«2n)[l/X  +  h/ (1-6)  -  T]/h 

+  «3  ♦  A  +  [o«'Ah/(l-.'Ah)l(V0D0  +  «j) 

+  VpDjj  +  «4iwu-e)  -  tj| 

♦  [i/x  ♦  h/u-M  -t  +  «d0*“x*7 

<l-*‘xh)  +  Dx]  (I) 

For  a  detailed  development  of  this  model*  see  H»t- 
gomery  [14).  Rote  that  given  a1#  a2,  *3,  a],  A,  a** 

1,  Vg,  Da*  D1#  and  tba  magnitude  of  the  proeass  shift, 
($)  may  be  easily  optimised  for  n,  h,  and  k  using 
direct  search  methods. 

Major  Applications 

Suppose  that  tha  quality  characteristic  of  in¬ 
tax  sat  la  represented  by  the  parameter  9,  and  that  x 
represents  the  soaple  statistic  corresponding  to  9. 
There  are  two  possible  values  for  9;  9  •  Bq  corres¬ 
ponding  to  the  in-eontrol  state  and  9  •  92  which 
represents  the  out -of -control  state.  Then  the  pro¬ 
bability  of  a  false  elans  la 

a  •  1  -  FCjLCL  £  x  S  VCl  I  •  *  9q)  (9) 

mi  tha  power  of  the  control  chart  la 

1-$«1-  FOCL  %  %  %  VCl  |  9  •  9t)  (10) 

whore  ICt  and  VCl  and  tbg  lower  and  upper  control 
limits  an  the  fhewhert  control  chart.  This  by  spe¬ 
cifying  the  quality  characteristic,  tha  sample  sta- 
f 1st is,  ami  tba  relevant  aampllng  distribution*  tha 
cost  nod  ala  of  tha  previous  section  could  he  applied 


to  any  type  of  Shewhart  control  chart.  For  example, 
if  the  quality  characteristic  of  interest  la  a  mea¬ 
surement  described  by  the  process  mean,  then  the 
x-chart  would  be  used;  thus  equations  (9)  and  (10) 
become 

a  *  29(-k) 
and 

1  -  8  «  +  *(-«yi-k) 

respectively*  where  9(z)  denotes  the  standard  normal 
distribution  function. 

The  major  applications  of  cost  models  such  as 
those  in  the  previous  section  have  been  to  the  x- 
chart,  the  x/R  chart  combination,  and  the  fraction 
defective  or  p-chart.  Some  work  has  also  focused  on 
procedures  for  the  simultaneous  control  of  several 
related  quality  characteristics,  and  the  treatment 
of  non-Shevrhart  control  charts  such  as  the  cumulative 
sum  control  chart  and  the  use  of  warning  limits  on 
the  x- chart .  Almost  all  reported  studies  assume  that 
the  process  failure  mechanism  is  exponential.  As 
noted  earlier*  this  is  a  critical  and  sometimes  un¬ 
warranted  assumption.  For  an  extensive  review  of  the 
literature  in  this  ares*  see  Montgomery  [14]. 

Most  of  ths  papsrs  cited  In  fl4]  contain  numer¬ 
ical  studies  end  usually ,  e  sensitivity  analysis. 

From  these  results,  it  is  possible  to  drew  certain 
general  conclusions  about  control  chart  design: 

1.  The  optimum  a  amp  la  site  is  largely  determined  by 
the  else  of  the  shift*  with  smeller  shifts  re¬ 
quiring  larger  samples.  Smell  shifts  say  require 
vary  large  samples*  possibly  as  larga  as  a*40  or 
more  on  «a  x  chert. 

1,  Chang ss  in  a^  and  a2»  the  fixed  and  variable  coata 
of  sampling,  affect  all  three  design  parameters. 
Increasing  the  fixed  cost  a.  increases  the  Inter¬ 
val  between  samples  and  leans  usually  to  slightly 
larger  samp lee*  while  increasing  the  variable  coet 
aj  usually  results  in  smell,  infrequent  ssaplss* 
bit  narrow  control  limits. 

3.  Ths  penalty  cost  for  out -of -control  production 
•1  mainly  Influences  ths  interval  between  sam¬ 
ples  h.  Larger  values  of  eg  Imply  emsllor  opti¬ 
mum  value#  of  h,.  A  similar  affect  is  observed 
by  increasing  X*  the  mean  number  of  proceee 
shifts  per  hour* 

4.  The  search  costs  S3  and  a«  mainly  effect  the  coo- 

!  trol  limits,  and  have  a  slight  affect  on  sample 

else.  Larger  values  of  S3  and  S3  result  in  wi¬ 
der  control  limits  end  lerger  sample  sixes,  as 
increasing  seer eh  costs  imply  that  fewer  false 
elsrms  ere  deelrable. 

9.  The  optimum  economic  control  chart  design  it 
relatively  insensitive  to  estimates  of  the  coet 
parameters.  The  coet  surfaces  are  usually  flat 
naar  the  optimum,  end  are  usually  steeper  near 
ths  origin,  so  that  it  is  batter  to  overestimate 
the  control  chart  design  parameters  than  to  un¬ 
derestimate  them.  The  optimum  design  is  most 


* 
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sensitive  to  errors  in  ost treat the  in- 
control  and  ouc-of -control  states  (or  the  mag¬ 
nitude  of  tht*  process  shift). 


Once  the  optical  n  and  k  are  determined,  the  Inter¬ 
val  between  samples  could  be  obtained  from 


This  last  finding  has  important  practical  impli¬ 
cations,  for  it  is  often  difficult  to  precisely  esti¬ 
mate  costs.  However,  the  mean  time  between  shifts 
and  the  magnitude  of  the  process  shifts  are  usually 
mora  easily  determined  from  process  performance 
data  or  from  the  engineer fs  knowledge  of  the  opera¬ 
ting  environment.  Note  that  one  should  exercise 
caution  in  using  arbitrarily  designed  control  charts 
and  empirical  "rules  of  thumb",  such  as  n* 5,  k*3  and 
h*l  for  the  x-chart.  In  some  cases,  particularly 
those  with  small  shifts  and  large  penalty  costs 
for  production  in  the  out-of-control  state,  signifi¬ 
cant  economic  penalties  may  result. 

A  Semi- Economic  Approach 

Instead  of  the  fully-economlc  modeling  approach 
outlined  above,  it  may  be  desirable  in  some  cases  to 
use  a  semi-economic  scheme  that  blenda  both  economic 
and  statistical  criteria  directly.  This  could  be 
useful  in  situations  where  the  analyst  la  unwilling 
or  unable  to  estimate  all  of  the  fully-economic 
mqd«l*«  cost  parameters,  or  where  he  wishes  a  com¬ 
promise  between  the  fully-economlc  and  purely  statis¬ 
tical  approaches. 

One  such  semi-economic  scheme  would  be  to  find 
the  values  of  n,  k.,  and  h  that  minimize  the  expected 
sampling  costs  per  cycle,  but  that  also  have  desir¬ 
able  average  run  lengths  for  the  control  chart  in 
both  the  ln-control  and  out-of-control  states.  Con¬ 
sequently,  we  would  Ilka  to  choose  n,  k,  and  h  so  as 
to 

minimize  t  •  (a^+a^n)  E(T)/h  (11) 

maximize  ARLQ  •  1/a  0.2} 

t 

minimise  ARI^  -  1/0-0)  03) 


h  * 


»  1 
aa3  +  ai  +  a^n  J 

Xa4(a-8)_1  -  0.5]  > 


(15) 


This  i9  an  approximation  for  the  optimum  sampling 
interval  h  given  n  and  k  suggested  by  Duncan  (7) 
and  Chiu  and  Wetherill  [A], 

These  formulations  of  the  optimal  control  chart 
design  problem  have  not  been  extensively  investi¬ 
gated.  It  would  be  of  interest  to  discover  hov  they 
compare  with  the  fully-economic  solutions,  and  with 
the  standard  designs  often  suggested  in  the  litera¬ 
ture. 


ECONOMIC  ASPECTS  OF  ACCEPTANCE  SAMPLING 


There  are  various  types  of  acceptance  sampling 
schemes;  single  sampling,  double  sampling,  sequen¬ 
tial  sampling,  continuous  sampling,  and  so  forth. 

The  quality  characteristic  being  Inspected  may  be 
either  an  attribute  or  a  variable.  For  various 
reasons,  more  attention  has  been  given  to  attributes 
sampling,  and,  in  particular,  single  sampling  plans 
for  attributes  have  been  studied  extensively.  We 
now  give  some  of  the  key  results  pertaining  to  the 
economic  design  of  these  plans  and  briefly  survey 
some  of  the  other  related  work. 

An  attribute  single  sampling  plan  is  indexed  by 
three  numbers:  the  lot  size  K,  the  sample  size  n, 
and  the  acceptance  number  c.  If,  in  a  random  sam¬ 
ple  of  site  n,  c  or  fever  defectives  are  found,  the 
lot  is  accepted,  while  If  more  than  c  defectives  are 
found  the  lot  is  rejected.  The  probability  of  lot 
acceptance  is 

c  /Nft\/N-Ne\  //N\ 

PC9)  •  l  (  i  ) (  n-d  )f  (n  )  (16) 


•here  ARX4  AHL^  are  the  average  run  lengths  in 
the  in-control  and  out -of -control  states,  respec¬ 
tively,  end  E(T)'  Is  glvsn  by  equation  (2) .  Recall 
that  a  and  1-0  sre  functions  of  the  sample  size  n, 
tbs  control  limit  factor- k,  and  the  magnitude  of 
the  shift.  Rota  that  this  is  a  mult  ler  Iter  Ion  opti¬ 
misation  problem,  and  would  generally  he  harder  to 
solve  than  the  optimization  problems  usually  asso¬ 
ciated  with  Che  economic  design  of  control  charts. 
However,  it  avoids  explicit  estimation  of  the  search 
coats  S3  and  e^,  and  the  penalty  cost  for  production 
la  the  out-of -control  state  04,  which  many  analyst t 
find  difficult. 

Another  possible  formulation  la  to  minimize  the 
stapling  coat  on  a  per  sample  basis  while  simulta¬ 
neously  mss  tailing  aRLq  end  minimising  ARLt .  Thus, 
<11)  becomes 

minimise  t"  •  a^  ♦  e^n  (1A) 

for  nbl,  sad  <12}  end  (13)  ere  unchanged.  This 
formulation  of  the  problem  Involves  only  n  and  k. 


where  6  le  the  lot  fraction  defective.  The  plot  of 
?(6)  versus  6  for  0  5  ft  5  1  is  called  the  operating 
characteristic  curve  (or  OC  curve)  of  the  plan. 
Traditionally,  we  design  acceptance  sampling  plans 
(l.e. ,  choosa  n  end  c)  so  that  P(ft)  peases  through 
or  near  two  points,  such  as  the  familiar  producer's 
and  consumer's  risk  points  (AQL,a)  and  (LTPD,B) . 

Whan  N  Is  large  relative  to  n,  the  hyper geometric 
distribution  in  (16)  esn  be  repieced  by  the  binomial. 

Prior  Distributions  for  Attribute  Sampling 

Lot  fraction  defective  8  Is  s  function  of  two 
sources  of  variability;  the  variability  of  p  the 
process  fraction  defective,  and  the  variability  of 
6  about  p.  It  la  usually  convenient  to  assume  that 
lot  quality  has  a  mixed  binomial  distribution;  that 
la,  each  lot  Is  produced  by  a  production  process 
that  ta  in-control  at  the  level  p,  but  p  varies  from 
lot  to  lot  according  to  a  probability  distribution 
f(p).  The  distribution  f (p)  Is  oftsn  celled  the 
prior  distribution  for  p  cr  the  process  curve.  It 
la  axtreaely  important  to  note  chat  If  the  process 
quality  Is  stable  such  that  £(p)  •  1  when  p  •  p(j 


and  f(p)  •  0  when  p  ^  pq,  th«n  *"here  is  no  need  for 
sampling  (this  result  is  called  Mood ' s  theorem;  see 

*>od  1153). 


where  a  is  d  constant  proportional  to  the  variable 
cost  of  sampling  and 


While  there  are  many  possible  choices  tor  the 
prior  distribution  for  p,  some  of  the  more  important 
are  the  following:  the  continuous  beta  distribution. 


?<p>  *  I  (d)pda-p)n"d 
d-0  v  ' 

is  the  probability  of  accepting  a  lot  of  quality  p.  ’  - 
Minimizing  L  with  respect  to  n  and  c  will  produce 
the  optimal  sampling  plan.  This  is  sometimes  celled 
the  Bayesian  approach  to  designing  a  sampling  plan. 


the  discrete  two-point  binomial, 

PCp-p^)  *  fi»  (18) 

and  the  normal-generated  distribution 

f(p)  -  cT1  explain2  -  *5(u-m)2/o2]  (19) 

where  p  *  ♦(-u) .  In  practice,  it  Is  important  to 
know  how  accurately  the  prior  distribution  must  be 
epee if led.  Generally,  the  analyst  does  not  possess 
sufficient  information  about  the  process  to  specify 
the  nr lor  with  great  confidence.  Fortunately,  most 
results  indicate  that  precise  specification  of  the 
prior  la  not  critical,  provided  that  a  reasonable 
distribution  la  chosen.  Continuous  prior  distribu¬ 
tions  are  generally  thought  to  be  more  appropriate 
than  discrete  ones,  end  the  beta  distribution  (17) 
has  been  used  extensively.  However,  when  the  under¬ 
lying  quality  characteristic  le  a  continuous  variable 
that  la  normally  distributed  within  each  lot  and  the 
mean*  of  this  quality  characteristic  also  has  a  nor¬ 
mal  prior,  then  the  prior  distribution  for  p  has 
the  form  (19) .  The  beta  and  normal -genera tad  dis¬ 
tributions  can  have  very  different  shapee  for  the 
same  mean  and  variance,  and  so  significant  differ¬ 
ences  In  the  optimal  sampling  plans  may  result.  For 
further  discussion  of  prior  distributions,  see 
Wether 111  and  Chiu  [21]. 


The  Major  A 
Attributes 


proaches  to  Sinils  Sampling 


Given  a  suitable  prior  distribution  and  e  set 
of  costs  or  losses  associated  with  sampling  plan 
operation,  It  la  desirable  to  choose  the  sampling 
plan  partmeters  that  minimise  the  total  cost.  Per¬ 
haps  the  most  widely  used  and  detailed  model  is 
that  of  Guthrie  and  Johns  [9],  A  simplified  form 
of  this  model  is  also  presented  by  Raid  [11].  The 
model  la  a  linear  cost  model.  All  linear  cost  for¬ 
mulations  lead  to  the  same  expected  loss  function. 
Held  utilises  the  concept  of  break-even  quality  pr , 
a  fraction  defective  value  at  which  It  le  just  as 
eoatly  to  accept  as  to  reject  the  lot.  The  expected 
lose  per  lot  Is 


lpt*J>)  f  1— »Cpl  ]t  lp)dp 

1  ,  V 

♦/  (p-Pr>Kp>f(pHp? 


The  major  focus  of  the  research  in  this  area  has 
started  with  the  Guthrie  and  Johns  model.  Hald  has 
been  a  major  contributor  in  the  field,  along  with 
some  of  his  coworkers.  The  major  emphasis  has  been 
on  finding  asymptotic  relationships  between  n  and  c 
for  various  process  curves,  and  in  producing  tables 
suitable  for  use  by  professional  practitioners.  Hald 
[10]  gives  a  number  of  significant  findings  in  his 
1960  paper.  One  part  of  the  paper  Investigates  the 
compound  hypergeometric  distribution;  that  is,  the 
probability  distribution  of  the  number  of  defectives 
d  in  a  random  sample  of  site  n  given  a  prior  distri¬ 
bution.  The  second  part  of  the  paper  is  directed 
towards  actually  finding  optimum  sampling  plans  for 
rectangular,  beta,  and  double  binomial  priors.  A 
general  solution  is  given  assuming  the  linear  cost 
model  and  inequalities  are  given  for  n  and  c.  In 
1965,  Hald  [111  provided  tables  for  the  double  bino¬ 
mial  prior  and  In  1968,  Hald  [12]  provided  tablea 
fox  the  beta  prior,  along  with  asymptotic  relation¬ 
ships  between  n  and  c,  and  between  N  and  n. 

Despite  the  significance  of  Held’s  work,  the 
tables  he  has  provided  are  often  difficult  to  use 
because  of  the  large  amount  of  information  required. 
An  alternative  approach  consists  of  formulating  an 
appropriate  cost  model  and  optimising  it  for  a  spe¬ 
cific  problem  using  direct  search  methods.  Consi¬ 
derable  work  in  this  area  has  been  done  by  G.K.  Ben¬ 
nett,  K.B.  Case,  and  J.W.  Schmidt  and  their  student*. 
For  example,  their  1972  and  1975  papers  [3] [17]  deve¬ 
lop  economic  models  for  single  sampling  plans  for 
dealing  simultaneously  with  multiple  attributes.  The 
cost  modsls  consist  of  a  component  representing  in¬ 
spection  costs,  a  component  representing  the  expec¬ 
ted  cost  of  lot  rejection,  and  a  component  represen¬ 
ting  the  Expected  coat  of  lot  acceptance.  Pattern 
search  la  used  for  model  optimisation.  In  general, 
direct  search  methods  are  a  very  effective  approach 
for  determining  economically  optimal  acceptance 
sampling  plant. 

The  economic  Impact  of  the  disposition  policy  for 
rejected  lots  hss  also  been  investigated.  While 
there  ere  a  number  of  possible  lot  disposition  pol¬ 
icies,  the  two  esses  that  have  been  Investigated  most 
extensively  ere  where* rejected  lots  are  either 
sera pp ad  or  screcnad  (100  percant  inspection).  For 
work  on  lot  disposition  policies,  sea  [18]  and [20]. 

Other  Work 

Thera  her*  baeu  many  other  studies  devoted  to 
the  economic  deelgn  pf  acceptance  sampling  plans; 
for  Instance,  a  1975  survey  paper  by  Wat her 11 l  and 
Chiu  [21]  cites  25)  references.  While  most  of  the 
work  focuses  on  single  sampling  for  at tributes,  there 


has  been  eexae  research  on  the  economic  design  of 
•ore  sophisticated  sampling  plans*  Foe  example, 
Stewart,  Montgomery  and  He ikes  [ 20 j  describe  a  proce¬ 
dure  for  the  selection  of  double  sampling  plans  for 
attributes  based  on  prior  distributions  and  costs. 
Models  are  presented  for  the  cases  where  rejected 
lots  are  either  screened  or  scrapped.  They  note  chat 
there  is  often  little  difference  in  cost  between 
economically  optimal  single  and  double-sampling  plans. 
However,  when  sampling  cost 3  are  large,  double  sam¬ 
pling  plans  have  much  to  offer.  They  also  observe 
that  arbitrary  double- sampling  plans,  such  as  those 
in  MIL  STD  105D,  may  be  very  far  from  economically 
optimum. 

The  effect  of  inspection  error  on  sampling 
plan  design  has  also  received  considerable  attention. 
If  an  Inspector  misclasslf let  good  and  bad  items  with 
constant  probabilities,  the  effect  is  to  translate 
the  0C  curve  of  the  sampling  plan,  so  that  the  actual 
or  effective  0C  curve  is  somewhat  different  from  the 
nominal  or  advertised  0C  curve.  If  the  probab lilt lee 
of  mlsclaeslf ication  are  known,  then  one  may  directly 
incorporate  thla  information  into  the  economic  design 
of  the  sampling  plan.  Generally,  the  presence  of 
inspection  errors  implies  that  larger  samples  are 
necessary.  A  good  review  of  the  literature  in  this 
area  Is  in  Dorris  and  Foote  16], 

Vary  littls  attention  has  been  given  to  accep¬ 
tance  sampling  by  variables.  Variables  sampling  is 
not  as  widely  used  In  practice  as  attributes  sampling 
for  several  reasons: 

1#  Variables  measurement  is  often  more  expensive 
(difficult)  then  attributes  measurement. 

2.  A  separate  plan  must  be  used  for  each  quality 
characteristic. 

3.  The  estimation  of  fraction  defective  assumes  s 
normally  distributed  quality  characteristic.  If 

.  this  assumption  is  violated,  the  tail  areas  may 

!  be  dramatically  aff acted  and  the  resulting  frac¬ 

tion  defective  estimate  grossly  in  error. 

4.  It  Is  possible  to  reject  a  lot  without  Actually 
finding  any  defectives,  and  this  often  upsets 
both  producers  and  consumer*. 

However,  remember  that  variables  sampling  can  greatly 
reduce  the  required  sample  else,  and  that  it  docs 
generally  provide  better  information  about  the  lot 
or  process  quality.  For  work  on  the  economic  design 
of  variables  sampling  plan,  see  Allor,  Schmidt  end 
Bennett  [l],  Scheldt,  Bennett  and  Case  118],  and 
Sclm&ldt,  Case  end  Bennett  [19].  Reference  [1]  deals 
with  the  situation  where  the  quality  character 1st ice 
•re  a  mixture  of  attributes  end  variables.  In  all 
of  these  studies,  the  approach  taken  is  to  formulate 
•  cost  modal  and  optimise  it  via  direct  search 
methods. 


The  last  20  years  have  seen  the  development  of 
numerous  techniques  for  the  design  of  process  control 
and  acceptance  sampling  schemes  based  on  economic 
considerations.  However ,  the  indication  is  that  very 


very  few  practitioners  have  Implemented  any  of  these 
techniques  (for  example,  see  the  surveys  by  Sanlga 
and  Shirland  116]  and  Chiu  and  Wetherlll  15]).  This 
is  surprising,  as  most  quality  assurance  managers 
claim  that  cost  reduction  and  increased  productivity 
Is  s  major  objective  of 'their  function.  In  many 
cases,  an  experienced  engineer  could  design  sn  appro¬ 
priate  technique,  perhaps  even  one  that  is  nearly  • 
economically  optimal,  but  the  use  of  a  formal  econo¬ 
mic  model  to  assist  the  analyst  Is  a  much  more  pre¬ 
cise  approach,  leaving  lees  to  judgment.  Further¬ 
more,  there  is  often  a  significant  economic  penalty 
associated  with  the  "standard",  judgment  designs  so 
frequently  used  in  practice. 

The  implementation  of  these  techniques  requires 
a  computer  program  of  the  cost  model  and  the  optimi¬ 
zation  procedure.  The  lack  of  availability  of  suit¬ 
able  computer  software  haa  certainly  slowed  practical 
implementation.  This  is  an  important  gap  between 
theory  and  practice  that  must  be  filled  before  the 
economic  design  of  quality  aeeursnee  techniques  will 
become  widespread. 

Many  practitioners  are  reluctant  to  use  these 
techniques  because  of  the  difficulty  In  estimating 
costs,  prior  distributions,  and  other  model  parame¬ 
ters.  Fortunately,  costs  and  prior  distributions  do 
not  heva  to  be  estimated  with  high  precision,  although 
some  other  parameters,  such  as  the  magnitude  of  the 
process  shift,  require  more  careful  determination. 
Sensitivity  analysis  of  the  specific  model  could 
help  the  analyst  discover  which  parameters  are  cri¬ 
tical  in  Ms  specific  application.  The  availability 
of  efficient,  interactive  computer  software  would  be 
of  elgcif leant  value  in  this  respect. 

.  We  will  Illustrate  the  use  of  a  simple,  inter¬ 
active  computer  program  for  the  optimal  economic 
design  of  sn  x-chart.  The  program  assumes  that  the 
Duncan  single-cause  process  model  is  appropriate, 
and  requests  the  user  to  Input  valuee  of  a^,  #2*  S3, 
*3*  *4 t  1*  6*  g,  end  D.  The  optimal  control  limit 
width  k,  interval  between  samples,  and  minimum  cost 
per  unit  time,  are  calculated  and  displayed  for  a 
range  of  sample  sizes,  along  with  the  n-risk  and 
power  of  tne  control  chart  for  each  (n,k,h)  combina¬ 
tion.  The  economically  optimal  control  chart  design 
may  be  found  by  inspection  of  the  cost  function  val¬ 
ues  to  find  the  minimum.  The  output  provided  enables 
the  Analyst  to  determine  the  sensitivity  of  the  coat 
surface  in  the  vicinity  of  the  optimum.  The  program 
Is  written  in  FORTRAN  for  a  CDC  CYBER- 7 4  computer, 
and  involves  less  than  100  lines  of  code.  Execution 
times  are  typically  one  second  or  lees.  Furthsr 
details  of  the  program  are  available  from  the  author. 

Consider  s  manufacturer  of  non-returnable 
glass  bottles  for  packaging  a  carbonated  soft  drink 
beverage.  The  wall  thickness  of  the  bottles  lr  an 
important  quality  characteristic.  If  the  wall  1* 
too  thin,  internal  pressure  generated  during  filling 
will  esuse  the  bottle  to  burst.  The  manufacturer 
has  been  using  x  and  R  charts  with  n*3,  k*3 ,  and 
taking  samples  every  h»2  hours  to  control  the  pro¬ 
cess.  He  wishes  to  compare  this  with  an  aconcmlcally- 
opt Intel  design  end  ss timet*  the  savings. 

Based  on  sn  analysis  of  the  quality  control 


technician**  salaries,  the  costs  of  the  test  equip¬ 
ment,  and  the  length  of  ti jze  requir'd,  it  in  eatt-**- 
ted  that  the  fixed  cost  of  taking  a  sample  is  Si. 00, 
and  the  variable  cost  of  sampling  is.  approximately 
$0.10  per  bottle.  It  takes  about  one  minute  (0.0167 
hours)  to  measure  and  record  the  wall  thickness  of  a 
bottle. 

The  process  is  subject  to  several  types  of 
assignable  causes  resulting  in  large  shifts,  typi¬ 
cally  of  about  two  standard  deviations.  (This  is 
why  the  "standard"  x-chart  design  n-5,  k*3  has  been 
used  -  it  is  known  to  be  reasonably  effective  for 
large  shifts).  The  previous  six  months  of  line  per¬ 
formance  data  Indicate  that  the  mean  time  between 
process  breakdowns  is  about  20  hours.  Thus  \»0.05  is 
the  parameter  of  the  exponential  distribution.  The 
average  time  to  investigate  an  action  signal  is  one 
hour.  Real  alarms  Incur  a  $25.00  search  cost,  on 
the  average,  vhlla  falsa  alarms,  which  are  more  dif¬ 
ficult  to  check  out,  incur  a  search  cost  of  $50.00. 

The  soft  drink  bottler  to  whom  the  bottles  are 
sold  has  a  policy  of  backcharging  the  bottle  manu¬ 
facturer  for  the  costs  of  cleanup  and  lost  production 
when  an  excessive  number  of  defective  bottles  burst 
during  filling.  The  past  six  months'  experience 
indicates  that  $100  per  hour  la  a  reasonable  estimate 
of  the  penalty  cost  of  operating  in  the  out-of -control 
state. 


Figure  1  snows  the  actual  computer  input  and  out¬ 
put  for  this  example.  Note  that  the  optimum  design 
has  n«5,  k*2.0y,  h*0.76  hours  (about  45  minutes),  and 
•  “inixnum  cost  of  $10.38  per  hn^r.  Thus  the  implica¬ 
tion  is  that  the  manuf^‘|-urer,s  existing  control  chart 
u’sign  is  good  with  respect  to  n  and  k,  but  that  he 
need**  to  samp,  more  frequently.  The  actual  cost  of 
his  cur^an#-  prc-edure  is  $12.05  per  hour,  or  about  ’a 
16  percent  penalty  c-jfit.  While  the  savings  per  hour 
seem  small,  this  is  a  continuous  production  process 
operating  three  shifts  per  day,  and  the  annual  sav¬ 
ings  are  over  $14,000.00. 

»  • 

After  looking  at  the  optimal  x-chart  design,  the 
bottle  manufacturer  suspects  that  he  may  have  incor¬ 
rectly  estimated  the  penalty  cost  of  out-of-control  . 
production  (a^),  and,  at  worst,  he  may  have  underes¬ 
timated  this  parameter  by  50  percent.  Therefore,  he 
reruns  the  program  with  a4«$150.00  to  investigate 
the  effect  of  mlsspeclfylng  this  parameter.  The  com¬ 
puter  output  f  shown  in  Figure  2,  indicates  that  the 
optimum  design  is  now  n-5,  V-2.99,  h-0.62,  with  a 
minimum  cost  of  $13.88  per  hour.  Note  that  the  pri¬ 
mary  effect  of  increasing  a$  by  50  percent  is  to  de¬ 
crease  the  interval  between  samples  from  45  minutes 
to  37  minutes.  The  program  could  be  used  in  this 
manner  to  quickly  and  easily  investigate  the  effects 
of  errors  in  specifying  any  of  the  model  parameters. 
Based  on  this  analysis,  the  manufacturer  elects  to 
adopt  a  45  minute  sampling  interval,  because  It  is 
administratively  simple  and  it  will  not  require  any 
additional  quality  control  technicians. 
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Figure  1.  Sample  Computer  Oitput 
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Figure  2.  Sample  Computer  <X»tput 
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